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Abstract 

■ In two-dimensional statistical models possessing a discretely holomorphic parafermion, we intro- 
duce a modified discrete Cauchy-Riemann equation on the boundary of the domain, and we show 

■ that the solution of this equation yields integrable boundary Boltzmann weights. This approach 
, is applied to (i) the square-lattice 0(n) loop model, where the exact locations of the special and 

■ ordinary transitions are recovered, and (ii) the Fateev-Zamolodchikov Zjv spin model, where a new 
^ , rotation-invariant, integrable boundary condition is discovered for generic N. 

oo 

(N ; 1 Introduction 

I Discretely holomorphic parafermions (also called lattice parafermions) are lattice analogs for the 

f holomorphic fields of Conformal Field Theory (CFT). They are lattice observables which satisfy 

a discrete version of the Cauchy-Riemann (CR) equation, and were recently identified in various 
two-dimensional critical lattice models, which are either spin models with a local interaction or loop 
models describing self-avoiding walks (SAWs) or extended interfaces. Generally, in spin models with 
Kramers- Wannier duality (see [T] for a recent review) , a lattice parafermion is constructed [U |3] 
by taking the product of lattice spin and disorder operators, providing an exact lattice analog for 
the operator product expansion (OPE) of the corresponding CFT [3]. In loop models, the lattice 
^ ' parafermion is expressed [SJ [5] in terms of a single open path attached to the boundary of the 

system. 

Lattice parafermions have found many applications to the rigorous study of lattice models 
(especially the Ising and SAW models) in the scaling limit, providing mathematical proofs for the 
exact results obtained by the Coulomb gas approach ^ i7j thirty years ago. In particular, the 
lattice parafermion is a crucial ingredient in the proofs [SJ E] of convergence of Ising domain walls 
I to the Schramm-Loewner Evolution (SLE), and in the calculation of Ising correlation functions [TO] . 

Also, for the SAW model on the honeycomb lattice, the lattice parafermion is essential to determine 
rigorously the location of the critical weights, both in the bulk [TT] and on the boundary [T^l [13]. 
This paper is concerned with lattice models at finite size, and addresses the relation between 

■ lattice parafermions and quantum integrability. Previous studies [21 El IS] have shown that many 
(spin or loop) statistical models admit a lattice parafermion exactly at their integrable point. In 
other words, the Boltzmann weights {M^(q;)} which satisfy the discrete CR equations on the rhombic 
lattice of angle a are also a solution of the Yang-Baxter (YB) equations, where a plays the role of 
the spectral parameter. In the present work, it is shown that a similar relation holds between the 
solutions of CR equations on the boundary and of Sklyanin's refiection equation [14] (or boundary 
YB equation). 

Two distinct models are considered in this paper. First, for the 0{n) loop model on the 
square lattice [IS] , following some ideas used in \TT^ , local CR equations around boundary faces are 
introduced, and their solution is related to the boundary integrable weights found in ^16) . When 
the rhombus angle is set to a = tt/3, one recovers the result of [H] for the honeycomb lattice. 
Then, we study the Fateev-Zamolodchikov Z^r spin model fl7j. We introduce a similar boundary 
CR equation, and show that its solution yields new non-trivial integrable weights, associated to 
free BCs. In both cases, we use the lattice parafermions constructed in [31[S] which satisfy bulk CR 
equations, and we define a new type of lattice domain where some of the boundary faces have a 
differnet shape from the bulk ones. 
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Figure 1: The local configurations of the square-lattice 0(n) model and the associated Boltzmann 
weights. 




(a) (b) 

Figure 2: (a) A domain of the rhombic lattice with angle a. (b) A domain Q of the rhombic lattice, 
where some boundary faces are replaced by triangles. In both cases, a is an arbitrary midedge of a 
boundary rhombus. 

2 The 0(n) model on the square lattice 

2.1 Integrable weights and Cauchy-Riemann equation in the bulk 

The 0(n) model on the square lattice was introduced in [TS]. It is a loop model, defined by the 
local configurations shown in Fig. [TJ and the Boltzmann weight of a loop configuration G: 

n[G] = X ^N^.iG) ^ 

where Ni{G) is the number of closed loops in G, and is the number of faces with configuration 
X. The partition function is denoted by Z . If the loop fugacity is set to 

n = ~2 cos4A , 

the integrable weights are given by 



t ~ sin(3 A — u) sin u + sin 2A sin 3A 

Ml = — sin(3A — u) sin2A 

U2 = — sinusin2A 

V — sin(3A — u) sinu 

wi — sin(3A — u) sin(2A — u) 

W2 — — sin(A — u) sin u , 



where u is the spectral parameter. We first consider a domain J7o of the rhombic lattice TZa (see 
Fig. [5^), and recall the results on the lattice parafermion for the 0(n) model. For a fixed edge a on 
the boundary, the parafermionic observable Fs is defined on the midedges {z} of the rhombi as [S] 

Fs{z):^^ E n[G]e-«(«), (2) 

Ger(a^z) 

where T{a — z) is the set of loop configurations comprising an open path between a and z, 0(G) is 
the winding angle of the path included in G, and s is called the conformal spin. For a polygon P 
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Figure 3: The two possible configurations of a triangular face contributing to ([8]) for a fixed outer 
configuration. 



with vertices given by the complex numbers (pi, . . . ,Pm) ordered counter-clockwise, we define the 
discrete contour integral of a function / around P as 

Y^mSz -p,)/ [^^^^^) , (3) 

where Pm+i '■— Pi- We define the discrete CR equations as 

Y,Fs{z)Sz^O, (4) 

o 

where the left-hand side stands for the discrete integral around any given rhombus of f^o- Fixing 
the outer loop configuration, this gives the linear system: 

t + ^Ul — ^T^^U2 — V = 

nui — TU2 — fJ.T^'^v + fi{nwi + W2) = 
— /ir^^Mi + fiTU2 + nv — T^'^wi — — , 



where we have set 



In [S] it was shown that the system ([S]) has a nonzero solution iff cos47rs = cosl2A, and one can 
choose the value 

. = ^-1. (6) 

TT 

Moreover, the solution of for u ^ is given by the integrable weights ([T]), if one sets the spectral 
parameter to 

u — {s + l)a = — a . (7) 

TT 

2.2 Boundary Cauchy-Riemann equation 

We now consider a domain fl where the interior faces are rhombi with a lower angle a, and the 
boundary faces are either rhombi or triangles (see Fig. [2]d). Moreover, we assume that the initial 
orientation of the open path at the boundary point a is horizontal. The Boltzmann weights on a 
triangular face are y,r (see Fig.[T]). 

Let us introduce the boundary CR equation: 



Re 



^F,(^)<5. 



(8) 



where the discrete contour integral is around a triangular boundary face T. The configurations 
G in ([2) which contribute to ([5]) are those where the open path ends on an edge of the boundary 
face T. Like in the bulk, one can fix the loop configuration outside T and sum over the internal 
configurations (see Fig. [3]). This yields 



Re 



,|[(s+l)a+7r(l-s)]^ _ [-(s+l)Q+7r(l-s)] J ^g^ 



We now have to choose an appropriate parameterisation for the spectral parameter. For this, 
we consider a particular instance of a domain fl: the infinite strip in a diagonal direction of TZa- 
In [16] . it is shown how to construct an integrable model on a strip with arbitrary, homogeneous 
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Figure 4: Left: the original square lattice L for the Zjv spins (full dots, full lines) and its dual 
lattice (empty dots, dotted lines). Right: the covering lattice TZ. 



spectral parameter in the bulk. If we denote i? and K the solutions to the bulk and boundary 
YB equations, the resulting model has bulk weights given by R{u), and boundary weights given by 
K{u/2). This suggests that in (|9]), one should set 

u= i(s + l)a 

instead of ([7]). In contrast, the conformal spin belongs to the definition of the observable Fg, and 
we keep its value The solution of ^ then reads 

r{u) = — cos^(3A + 2w) 

(10) 

y{u) = cos ^(3A — 2-u) . 

These are exactly the boundary weights identifed in |16| . which satisfy the boundary YB equation 

Ri2{u - v)K2{u)Ri2{u + v)K2{v) = K2{v)Ri2{u + v)K2{u)Ri2{u - v) . 

In [TB], it is argued that PH)) corresponds to the special surface transition in a certain regime of A. 
Under the change A — > A + tt, the bulk weights are invariant (up to irrelevant sign factors), whereas 
the boundary weights correspond to the ordinary surface transition. 

Note that, if we specialise to m = A (a = ^), we recover the critical bulk and boundary weights 
of the 0(n) model on the honeycomb lattice, which were found through a type of boundary CR 
equation in [T^ for a specific domain $7. 

3 The Fateev-Zamolodchikov model 



3.1 Definitions 



The Zat model consists of spins aj living on the sites of the square lattice C (full dots in Fig. H]), 
and taking the values: 

, AT" I, / ^ZTT \ 

o-j G {l,w,...,cj }, where w = exp ( — I . (11) 
We shall then write aj = oj^^ , with Vj G Z. We also introduce the crossing parameter 

2N 



The Boltzmann weight of a spin configuration is: 



rrr 1 TTw/ \ u Tuif \ j^i'^l'^i - "^j) if (ij) is horizontal, 

n[a] =||VF(a„(7,), where Wia^.a^) ^ i ' ^ -f/ A- i 

■^^ I W[X — u\ri — Tj) it [ij) IS vertical. 



(12) 



where u is the spectral parameter. Moreover, we restrict to a periodic and symmetric function 
W{u\r): 

Wiu\r + N)^W{u\-r)^Wiu\r). (13) 



The partition function is then defined as 



Z = J2^[a], (14) 
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and a typical correlation function of n spin operators is of the form 

/ ki ko kn\ ^ \ ^ Trr 1 1^1 '^^ 



where {ji, . . . , jn} are n sites of the lattice and . . . , mn\ are n arbitrary integers. As a 
consequence of the symmetries \Vi\ of the Boltzmann weights, the model enjoys Kramers- Wannier 
duality, and in particular one may define dual spin operators, also called disorder operators /i'™-*. 
These live on the dual lattice (empty dots in Fig. |4|), and are defined through their correlation 
functions. For the purposes of the present study, it is sufficient to define correlation functions 
involving two operators. Let fc and i be sites of the dual lattice and an oriented path 
from fc to £ on C! . One then defines: 

(A.i-'"Vf^(-.-)):=^E n m^-'^.) n m^.^'"^.) (•■•), (15) 

{<t} {^])(i-ltl 

where the dots denote any function of the spins, 7^ is the set of edges of C which cross 7fe^, and 
we have taken the convention that ct^ (ctj) is on the left (right) of ^m- The symmetries p3p imply 

that is independent of the choice of the path ^m- 

In the following, we shall consider a deformation of the square lattice C into a rectangular lattice 
of aspect ratio tan ^- The union of C and its dual is then a rhombic lattice of angle a, called the 
covering lattice and denoted TZa ■ 

Later in this Section, we will use the discrete Fourier transform 

r=0 r=0 

3.2 Integrable weights and Cauchy-Riemann equation in the bulk 

The integrability condition for a periodic system {i.e. the commutation of transfer matrices on the 
cylinder) has the form of a star-triangle equation 

Af-l 

W{X ~ u\ri ~ r')W{u ^ v\r2 - r')W{\ - v\r-i - r') 

r'=0 

= W{u\r2 - r3)W{X ~u~ v\ri ~ ry,)W{v\ri - , (17) 

which can be formulated as the Yang-Baxter equation on the covering lattice (see Fig. [5]), and a 
solution to this equation was given in p/7j : 

W(u|0) = l, m^|r) = n '"lift mil forrG{l,...,iV-l}. (18) 

pJ-Q sm[(2p + 1)A + u] 

This solution is self-dual, in the sense that 

W{u\k) = W{X-u\k) . (19) 

The construction 3 of the discrete parafermion in the Zjv model is the lattice version of 
the parafermionic current in the Zjv CFT 4 . In the lattice model, as well as in the CFT, the 
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Figure 6: A face of the rhombic covering lattice TZa in the Z^r model. The dotted line represents an 
arbitray path 7 on C from /u(a) to /ii. 



parafermion ip is obtained as the product of a spin operator cr and a disorder operator fi on 
adjacent sites. Consider a point z on the middle of an edge of TZa, and denote j{z), k{z) the sites 
of C,C' adjacent to z. We shall use the short-hand notations 

a{z) := , fi{z) := , fi*{z) := . 

The lattice parafermion ^{z) is defined as |2| 

ij{z) a{z) . (20) 

The observable Fs{z) is then defined by fixing a point a on the boundary, and setting 



^ — is6{a.z) 



riaWz)) 
a*{a)cr{z)fi*{a)fi{z)) , 



(21) 



where 6{a, z) is the angle between the vectors k{a)j{a) and k{z)i{z), and s is the conformal spin. 
As in the 0(n) model, one defines the discrete CR equation on TZa as 



Y,Fs{z)5z = 0, 



(22) 



which can be written (see Fig. [SJ 

e ^ ' (cr (a)(Ti^ (a)^i [(e -e 0^02) + {e ^ -e cr^ 0-2)^1^^2] t) , (2d) 



where 9{a) is the angle between k{a)j{a) and the horizontal, and we have set 

u^^{l-s){TT^a). (24) 

From the definition of fi, reads 

i^a*(a)ai J] Wia.,a,) [] ^F(a„wa,) 

{<t} (ii)^7^U{(12)} (ii)67-L 

X [(e-*" - e™a*(T2)PF(ai,a2) + (e™+"^ - e-™-*"ViV2)VF(ai, was)] =0, (25) 

where 7 is an arbitrary path on C going from fi{a) to /ii (see Fig. IH]). The bracket in may be 
written as 

/(r) (e^™ - e*"w'')PF(M|r) + (e*"+*^^ - e-™-*'^'*w'-)W^(u|r + 1) , (26) 

where 

Cr*tT2 = w"" . 

Therefore, a sufficient condition for ([25)1 to hold is that /(r) vanishes for all r, which yields the 
recursion relation 



sm 



W{u\r + 1) = W{u\r) x (-e'''^+-) 



(2r+l)7i- 
2Af 

(2r+l)Tr 
2N 



(27) 
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Figure 7: The two-row transfer matrix along one direction of the covering lattice, with open boundary 
conditions. Thick lines denote interaction terms between the spins. 




Figure 8: The boundary Yang-Baxter equation for the Z^v model. Thick lines denote interaction 
between the spins, and spins joined by a double line must be equal. 



Compatibility of ((27|) with the symmetries (IT3| fixes the value of the spin 



^ = 1-^, (28) 

and for this value, the solution of ^7} is given by the integrable weights (fT5)) . 

As explained in j3], each of the (N — 1) discrete parafermions i{j'-"^\z) = (T™(z)/i^™)(z) with 
charge m G {1, . . . ,N — 1} satisfies the CR equations, but for a different function VF*^"'-', and a 
conformal spin s„i = m{N — m)/N. In fact, VF^™) is simply related to W by the global 
transformation aj cr™. Hence, one can restrict to the case m = 1 exposed above, without loss 
of generality. 



3.3 Boundary Yang-Baxter equation 

Following |18) . we consider the boundary conditions which ensure the commutation of two-row 
transfer matrices (see Fig. [T]). We denote yR,L(o'i, fj) the Boltzmann weight on the right and left 
boundaries, and the partition function on the strip reads 

Zstnp = J2 n ^i'^^^'^j) n YR{a,,a,) Yl YL{a,,aj). (29) 

{ct} (z_;)£bulk (zj) bright bound. (ij)f^lcft bound. 

In [TB], fixed-spin integrable boundary conditions {i.e. Y{ai,aj) ~ 6ri,a5rj,i3 with a, /3 fixed) were 
considered, and the associated surface free energies were obtained. In the present work, we address 
the boundary YB equation for rotation-invariant boundary conditions, i.e. an interaction of the 
form 

Y{cj,,<j,) = Y(u,£,\n~r,), (30) 

where ^ is a complex parameter associated to each boundary [boundary field), and Y satisfies the 
symmetries p^ . The boundary YB equation then reads (see Fig. [5]) 

W{u - v\ri - r2) ^ Y{u, ^\ri - r')W{u -h v\r2 - r')Y{v, ^jrs - r') 

r'=0 

= W{u - v\r2 - rs) ^ Y{v\ri - r')W{u + v\r2 - r')Y{u, ^[rg - r') , (31) 

r'=0 
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Figure 9: A boundary face in the Z^y model. Like in Fig. [6l the dotted hne represents an arbitray 
path 7 on £' from /i(a) to ^i. In the discrete contour integral ([33|) . only the edges represented by 
single lines contribute. 



or, in Fourier space: 

v\k + £ — m) 



N-1 

= Y{v,S,\e) ^ t?(u-z;|TO)vi7(u + i;|fc + £-TO)y(u,^|/c-TO). (32) 

These functional equations for Y are difficult to solve for general N, but we shall now show that a 
nontrivial solution can be found by use of the boundary CR equation. 

3.4 Boundary Cauchy-Riemann equation 

We consider boundary faces in the shape of a pentagon, determined by the angles a and /3 (see 
Fig. (HI . The contour integral of around such a face P reads: 

J2Fsiz)Sz = e-^''^'^) (^*(a)(TiMi [(g-^™ _ e^^a^Vs) + (e'«+"^ - e-^«-^^ViV2)M>2] ) , (33) 
p 

where we have selE] 

u = ^{1 - s)i7r - a) , and ^ = ^{1 - s){7r - p) . (34) 
Like for the bulk CR equation, using the notation crj['(T2 :— lj^ , p3p is rewritten as 

^F,(z),5z= ^^^a*(a)ai [] W{<j,,aj) [] VF(a„ wa,) x J(r) , 

P {a} fe)^7^ fe>e7^ 

where J(7') is defined as 

J(r) := (e-2™ - e2™a;'^)y(u, ^|r) + (e'«+''^^ - e-*«-''^''cj'^)F(u, + 1) , (35) 
and its Fourier transform is 

J{k) := (e-2» + e'^+''''uj-'')Y{u, ^\k) - (e^™ + e-'^-''"'uj-''-^)Y{u, £_\k + 1) . (36) 
In analogy with the 0(n) model, we introduce a boundary CR equation of the form 



Re 



e»¥'(fc)j(fc) =0, (37) 



where (p{k) is real. It turns out that Y only satisfies the symmetry relations (jl3p for the choice 
(p{k) — (4fc + 2)A, and the solution of ([57)1 then reads 

{1 for fc = 0. 

"y^ Sin [{2£ + 1)A - u + sin [(2£ +l)X-u-^] (38) 
ii sin[(2^+l)A + u + C]sin[(2£+l)A + u-C] ^ t^' • ' • ' ^b 

A remarkable fact is that the weights p8l) are a/so a solution of the boundary YB equation p2p . 
The solutions associated to the other parafermions ?/;(™) also satisfy the boundary YB 

equation, and they are related to (fT8|) - ((38l) by a simple transformation aj — ?> cr™. 



^Like for the 0(n) model, for a given angle a, the spectral parameter for boundary weights is half the bulk value (|24|l . 
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3.5 Physical interpretation of the integrable boundary condition 



To indentify the BC corresponding to (j38l) . we consider the diagonal-to-diagonal transfer matrix, 
which is constructed from the two-row transfer matrix by the procedure described in (I6j . If the two- 
row transfer matrix with horizontal spectral parameter u, vertical spectral parameters (vi, . . . , vn) 
and boundary parameters ^l,r is denoted 

t{u,^L,£,R]vi, ■.■,Vn), 

then the diagonal-to-diagonal transfer matrix is given by 

/ u 



u u u u 



In this setting, the bulk weights are given by for horizontal edges and W{X — u\-) for vertical 

edges, whereas the left and right boundary weights read 

Y{\-u/2,iL\r)Wi\-u\r), and Yiu/2,^R\r). 

The latter coincide for the value — £,r — u/2, where both are proportional to W{X — u\r). Hence, 
for this choice of boundary parameters £,l,£.r, the boundary weights are equal to the bulk weights 
in the vertical direction, and thus psp corresponds to free boundary conditions for the spins {(Jj}. 



4 Conclusions 

On the two examples we have studied, we have shown that the integrable boundary weights can be 
obtained by imposing a simple linear condition - which we have called the boundary CR equation 
- on the discrete contour integral 

p 

where P stands for a modified boundary face. In the case of the 0(7i) model, we have recovered 
the integrable boundary weights which were obtained |16| through the mapping to the 19-vertex 
model, whereas for the Zjv model, we have found new integrable boundary weights. However, unlike 
the bulk CR equations, where the conformal spin s is generally obtained by simple consistency 
conditions, it is not clear yet how to extract conformal data from the boundary CR equations. 

Our approach is quite general, and is likely to extend to any lattice model where a solution 
of the bulk CR equation has been identified. Also, an interesting continuation of this work would 
be to study the boundary CR equations in loop models with "non-trivial" integrable boundary 
conditions, i.e. where loops touching the boundary get a different weight 19'. 

Finally, an important possible application of our results would be to set up proofs of convergence 
to variants of SLE as in [5112], in the presence of various integrable boundary conditions. 
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